Motifs haven’t been updated in this yet ref = “total_ratio” changed to “self” when using show_burden_for_selected_features Use a newly calculated set of motif change scores to define motif SNVs, but only focusing on motifs that more relevant Here motifs are defined as relevant if they are 1) in top 50% brain expressed genes, and 2) are enriched in at least one gain-enriched brain module Partition enhancers based on distance to TSS based on Yanyu’s pipeline. Basically, her pipeline takes intersections between enhancers and distance ranges, and remove bases of enhancers that are not overlap.
non ASD genes are the bottom 20% not plotting mutation is 5’UTR or 3’UTR added one row for showing the union of enhancers and promoters Only have Fantom, Noonan and Whole genome Look at more genesets.
getwd()
## [1] "/media/yuwen/F/ASD/analysis"
knitr::opts_chunk$set(warning=FALSE, message=FALSE)
## [1] "/media/yuwen/F/ASD/analysis"
Whole genome
Noonan enhancers
Fantom fetal brain
Brain roadmap adult brain
Brain roadmap + Noonan brain adult brain
plot_summary = plot_summary[!is.na(plot_summary$gene_list),]
plot_summary = plot_summary[!duplicated(plot_summary),]
#scaling factor for plotting, which is the burden of baseline. Used for 160229 data only
scaling = 0.73